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In systems displaying a bulk first-order transition the order parameter may vanish continuously 
at a free surface, a phenomenon which is called surface-induced disorder. In the presence of surface- 
induced disorder the correlation lengths, parallel and perpendicular to the surface, diverge at the 
bulk transition point. In this way the surface induces an anisotropic power-law singular behavior 
for some bulk quantities. For example, in a finite system of transverse linear size L, the response 
functions diverge as ]r j ( d - 1 ) z + 1 ^ w here d is the dimension of the system and z is the anisotropy 
exponent. We present a general scaling picture for this anisotropic discontinuity fixed point. Our 
phenomenological results are confronted with analytical and numerical calculations on the two- 
dimensional g-state Potts model in the large-g limit. The scaling results are demonstrated to apply 
also for the same model with a layered, Fibonacci-type modulation of the couplings for which the 
anisotropy exponent is a continuous function of the strength of the quasiperiodic perturbation. 



-a 

G 

O 

(N 
> 
OS 

m 
o 



i 

C 

o 
o 



X 



I. INTRODUCTION 

At a first-order phase transition several physical quan- 
tities, such as the order parameter and the energy density, 
have a discontinuous behavior in the bulk and the cor- 
relation lengths remain generally finite. In spite of the 
absence of a diverging "true" length scale in the problem, 
a successful scaling theory for first-order transitions has 
been developed, 1 ' 2 with scaling exponents such as 



(3d = 0, old = Id 



1. 



7D 



(1.1) 



which are compatible with discontinuities. The divergent 
length scale appearing in the scaling theory is associated 
with the persistence length of the coexisting phases, 2 be- 
having as |i|~' yD for a temperature (field) devi- 
ation from the transition point, with critical exponents 



I'D = 



l/d 



(1.2) 



where d denotes the dimension of the system. In the 
renormalization group analysis of specific models, first- 
order transitions have been identified, 1,3,4 which are con- 
trolled by a so-called discontinuity fixed point (DFP) 
with thermal and magnetic scaling dimensions yt — 
In At/ In b — d and yh — In Ah/ In b = d. These di- 
mensions are compatible with the value of vu and v D 
in Eq. (1.2) whereas other scaling exponents in Eq. (1.1) 
follow from scaling and hyperscaling relations. 

For a system with a free surface at a bulk first-order 
transition two diverging lengths can be identified, £|| as- 
sociated with the correlations in the direction parallel 
to the surface and £j_ in the perpendicular direction. 



Due to missing bonds, the surface is more weakly cou- 
pled than the bulk. As a consequence, the surface or- 
der parameter vanishes continuously at the bulk tran- 
sition point. This phenomenon, called surface-induced 
disorder 5 (SID), is closely related to wetting phenomena. 6 
It has been studied using mean-field theory 7,8 and Monte 
Carlo simulations 9 and has been observed in some exper- 
imental systems. 10 For the two-dimensional (2D) Potts 
model exact results have been obtained in the large- 
q limit and the universality of the critical and tricriti- 
cal exponents in the bulk first-order transition regime, 
q > 4, has been accurately demonstrated by density ma- 
trix renormalization calculations. 11 

The presence of diverging length scales in the problem 
suggests a possible scaling behavior for singular quanti- 
ties, both at the surface and in the bulk of the system. 
Here we notice two peculiarities concerning the bulk scal- 
ing relations. First, in the geometry with a free surface 
there is an inherent anisotropy in the problem since the 
two length scales are related by 



£il ~ £5 



(1.3) 



with an anisotropy exponent z which is generally greater 
than 1. Consequently, scaling for SID is strongly 
anisotropic in contrast to the original DFP theory, which 
was developed for isotropic scaling in the bulk. 2 Our sec- 
ond remark concerns the behavior of the system as the 
thermodynamic limit is approached. In a system of linear 
size L, the disordering effect of the surface in the bulk is 
limited to a temperature range At ~ L^ x l Vl - around the 
critical temperature, where v±_ is the perpendicular cor- 
relation length exponent [see below in Eq. (2.2)]. In the 
thermodynamic limit, when L — > oo, the SID is present 
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in the bulk just at the transition point, when the scal- 
ing functions are singular, having either a jump or 5- 
function-like singularities. In a finite system, however, 
the scaling functions behave regularly, so that one has a 
"true" finite-size scaling behavior, which is measurable, 
for the specific heat or the susceptibility at the transition 
point. 

Another problem of theoretical interest at a first-order 
transition point accompanied by SID is the influence of 
different types of perturbations, such as localized or ex- 
tended defects, 12 quenched disorder, 13 etc., on the bulk 
and surface transitions. Here we mainly consider the ef- 
fect of a quasiperiodic or, more generally, aperiodic mod- 
ulation of the couplings in the bulk. The influence of 
such perturbations on order-disorder phase transitions 
has been studied in detail for the Ising model, 14 in par- 
ticular in 2D with a layered structure (i.e., when the ape- 
riodic perturbation has a ID variation) and in the related 
quantum Ising spin chains. Among others, for marginal 
aperiodic sequences, coupling-dependent anisotropic crit- 
ical behavior has been found. 15-17 

At a first-order transition point, such as for the 2D 
Potts model with q > 4, 18 ' 19 the relevance or irrelevance 
of the effect of aperiodic perturbations is unclear. One 
can use an analogy with the rigorous result of Aizcn- 
man and Wehr, 20 stating that in 2D a first-order tran- 
sition always softens into a second-order one, due to 
quenched disorder. Such a softening of the transition 
was indeed observed in a Monte Carlo study of the 2D 
Potts model with q > 4 for some types of layered ape- 
riodic perturbations. 21 However, due to computational 
difficulties, quantitative results about the phase transi- 
tions, such as critical exponents, could not be obtained. 

In the present paper we use another approach to study 
the effect of aperiodic perturbations at a first-order tran- 
sition point and to obtain quantitative results. For this 
purpose we consider the 2D Potts model 19 on layered 
lattices and take the extreme anisotropic limit, 22 when 
the logarithm of the transfer matrix of the model cor- 
responds to the Hamiltonian of an aperiodic quantum 
Potts chain. Then we consider the large-g limit of the 
model, for which a consistent l/^/q expansion is defined 
and the leading term can be treated by analytical or very 
efficient numerical methods. In analogy with the results 
on the marginal quantum Ising chains, one expects that 
with marginal perturbations the critical exponents, in 
particular the anisotropy exponent, become coupling de- 
pendent. Indeed, with a Fibonacci-type modulation of 
the couplings, which is related to the 2D Penrose quasi- 
lattice, we have found continuously varying anisotropic 
scaling. In this way we have obtained a more general 
testing ground for the scaling relations in the presence of 
SID. 

The structure of the paper is the following. In Sec. II 
we present the phenomenological scaling theory for SID 
both in the bulk and at the surface. In Sec. Ill we in- 
troduce the quantum Potts model and its representation 
in the large-g limit. In Sec. IV solutions of the model on 



the regular lattice and on the Fibonacci quasilattice are 
confronted with predictions of the SID scaling picture. 
The results are discussed in Sec. V. 



II. SCALING BEHAVIOR WITH 
SURFACE-INDUCED DISORDER 

We consider a d-dimensional (magnetic) system in the 
slab geometry with two surfaces at a large finite distance 
L apart, whereas the extent of the system is infinite in 
the other (d — 1) dimensions, parallel to the surface. The 
system displays a bulk first-order transition and has q 
degenerate ordered phases. The state of the system on 
the surface at I = L is fixed into one of these phases 
(by applying a local magnetic field in the appropriate di- 
rection) whereas the other surface at I = 1 is free. The 
temperature t is measured relative to the transition point 
and a bulk ordering (magnetic) field of strength h is ap- 
plied. Due to the disordering effect of missing bonds, the 
order parameter profile mi is monotonically decreasing 
towards the free surface. The surface order parameter 
mi may behave in two qualitatively different ways as the 
temperature is increased to the transition point. Either 
mi , like the bulk order parameter m^ — m L / 2 , has a finite 
jump at t = (first-order surface transition) or — and this 
is the case considered below — mi vanishes continuously 
as 

mi~(-t) ft , t<0, (2.1) 

with a surface magnetization exponent /3i (second-order 
surface transition) . In this case the extent of the surface 
region where mi ~ mi (1 < I < is diverging at 
the transition point as 

a ~ (-ty v± . (2.2) 

As mentioned in the Introduction, this phenomenon as- 
sociated with the long-distance disordering effect of a free 
surface is called SID. 

In the region of SID, correlations between spins in the 
direction parallel to the free surface decay exponention- 
ally, which defines a new correlation length with an 
asymptotic temperature dependence 

£|| ~ (-*)-"" • (2-3) 

According to Eq. (1.3) the correlation-length exponents 
satisfy the relation v« = zv±. Going through the tran- 
sition point by varying h at t = leads to a similar 
behavior with —t replaced by h in Eqs. (2.2) and (2.3), 
while the corresponding critical exponents are v\ and vfr, 
respectively. 

Now using the scaling hypothesis for strongly anisotro- 
pic systems 23 we can write the following relation for the 
singular part of the bulk free-energy density: 

f(t, h, L) = b-W-^+^fitb 1 ^ , hb x ^,L/b) , (2.4) 
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when lengths are rescaled by a factor b > 1. By deriva- 
tion one obtains a similar scaling relation for the bulk 
order parameter, 

df 

m b (t,h,L) = — 

= b-^-^+^'^mbitb^M'^^/b) , (2.5) 
and for the excess internal energy per degree of freedom, 

df 



Au(t, h, L) = u(t, h, L) - u(0, 0, oo) 



dt 



= b-l(d-l)z+l]+l/v ±Au ( tb l/v ± ^ hb l/v'l jL / b } ( 2 g) 

Since the first-order transition is accompanied by a mag- 
netization jump and a finite latent heat, both mb(t, h, L) 
and Au(t, h, L) approach a scale-independent, finite limit 
when L — ► oo. Thus Eqs. (2.5) and (2.6) lead to the scal- 
ing relations 



1 



(d-l)z + l 



(2.7) 



One may notice that in the case of isotropic scaling, i.e., 
with z = 1, we recover the result of the DFP in Eq. (1.2). 

For the susceptibility \ — d 2 f/dh 2 and the specific 
heat c v = d 2 f /dt 2 , using the scaling relation in Eq. (2.4), 
we obtain 

X (t, h = 0,L) = t-^Xi(Lt^) = L ( - d ~ 1 ^ z+1 x 2 (Lt" ± ) 
c v {t, h = 0,L) = t- a °c vl (Lt^) = L^ z+1 c v2 (Lt^) , 

(2.8) 

with the DFP exponents jd and old in Eq. (1.1). The 
scaling functions behave as Xi(y) ~ c vl (y) ~ y( d ~ 1 ) z + 1 
and X2(y) ~ c v2 (y) — const as y — > 0. Thus, in a finite 
system, when approaching the transition point, one can 
observe the development of a singularity with the DFP 
exponents in Eq. (1.1). However, in the thermodynamic 
limit L — > oo, this singularity transforms into a 5 function 
as it should at a first-order transition point. Indeed the 
temperature range where the surface can influence the 
behavior of the bulk is limited to At <~ L~ 1 / u± , so that 
Atx(t = 0, h = 0, L) w 1, and in this way x{t, h = 0, L) — 
and similarly c„(t, h = 0, L) — is a representation of the S 
function as L — > oo. On the other hand, at the transition 
point the response functions show a power-law finite-size 
dependence, with the same critical exponent (d — l)z + 1 
as given in Eq. (2.8) and which should be measurable in 
finite samples. 

For the correlation lengths one can similarly write the 
scaling relations 

U(t,h,L) = biy_(tb x l v ^,hh x l v ^,Llb) , 

£H it, h, L) = V- £|| (tb 1 ^ , hb 1 ^ , Lib) , (2.9) 

from which the temperature (field) dependence in 
Eqs. (2.2) and (2.3) follows by taking the length scale 

aab = t- v± {b = h-<). 



We close this section with a similar analysis of the sin- 
gular behavior at the surface. For the singular part of 
the surface free energy per degree of freedom we write 
the scaling relation 



fi(t,h,hi,L) 
= b-^-V* hitb 1 '^ , hb 1 '^ , M (d ~ 1)2 



,L/b), (2.10) 



where a surface ordering field hi with anomalous scaling 
dimension x\ has been included. For the surface order 
parameter m\ = df\/dh\, Eq. (2.10) leads to the follow- 
ing scaling behavior: 

mi(L) ~ L~ Xl , mi (t)~t XlU± . (2.11) 

Comparing Eqs. (2.2) and (2.11) we obtain 



Pi = XlV±_ 



(2.12) 



Taking other derivatives of fi one obtains similar scaling 
relations for the surface susceptibility xi,i = d 2 fi/dh 2 , 
for the excess surface order parameter m s = dfi/dh, etc. 

To conclude this section, let us remark that at a first- 
order transition point which is accompanied by SID, the 
bulk singularities in finite systems involve the anisotropy 
exponent z and the dimension d of the system whereas 
the surface singularities involve a new exponent xi. In 
the following section our scaling assumptions are con- 
fronted with analytical and numerical results on specific 
models. 



III. THE TWO-DIMENSIONAL POTTS MODEL 
IN THE LARGE-q LIMIT 

We start with the Hamiltonian of the classical g-state 
Potts model on the square lattice, 19 

oo L 

-? H = E X^i(0[<*Kk+i-*j,fc)-i] 

k— — oo l—l 

oo L — l 



5(s l+ i tk - S; ifc ) - - 



(3.1) 



where si_k — 1, 2, . . . , q is a g-state Potts variable defined 
modulo q at site (I, k). The vertical [horizontal] couplings 
Ki(l) [K 2 (l)] are ferromagnetic and may depend on the 
horizontal position I. Constant terms have been added in 
order to simplify some expressions in the following. The 
transfer operator T has matrix elements 



<*|7V> = (s\U\s}{s\V\s') 



■ L-l 



5(si +1 - si) - - 



(s\U\s)=exp\J2 K ^ 
{ i=i 

(s\V\s')=e^^2Ki(l)[8(s' l -s l )-l]^ , (3.2) 
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between the Potts states |s) and \s') associated with two 
successive rows k and k + 1. Let us associate with each 
site / along a row the Potts operators f2; and M; such 
that 24 

Sll\si) = exp (i^-s^J \si) , fl\ = Clf 1 , 
Mt\ Sl ) = \si + 1) , M\ = M- X . (3.3) 

Thus fli is diagonal in the basis of the Potts states 
whereas M/ is a ladder operator. On the same site, they 
satisfy the commutation rules 

.2tt 



(3.4) 



Mifli = exp ( -i— ) QiMi 



M}n t = exp | /— | \I, 



2tt 

q 



t 

i ' 



and commute on different sites. Making use of the iden- 
tity 



1 9-1 

5(si +1 - si) = -J^cxp 



p=0 



. 2pw 

i (s/+i - s;) 

q 



(3.5) 



the diagonal operator U can be written as 

^ = exp{E^B^t +1 ) P \ • (3-6) 
I ;=i g P =i J 

In Eq. (3.2) the diagonal matrix element (s|F|s) is equal 
to 1. When the two states \s) and \s'} differ only on site 
I, we have (s|U|s') = exp[— K\{1)\. Thus one may write 
the nondiagonal part of the transfer operator as 

V= l + ^expH^OE Mf + ..., (3.7) 
i=i p=i 

where the ellipsis stands for higher-order terms involving 
products of exponentials coming from different sites. 

In the extreme anisotropic limit 22 where K\(l) — > oo, 
K2H) — > 0, and 1? — > in such a way that the ratios 



Ji 



expJ-^Q] 



(3.8) 



remain finite, the operators in Eqs. (3.6) and (3.7) are 
now given by 

u = 1 + $ E 7 E ( n ' n w) p + °^ 2 ) > 

;=i * P =i 



£ r 9-1 

y = 1 + 1?^ — J2K + °(d 2 ) ■ 
1=1 q P =i 



(3.9) 



To first order in U and V commute and the trans- 
fer operator can be rewritten as the evolution operator 



T = exp(— dTi) for the infinitesimal Euclidian time step 
■d. According to Eqs. (3.2) and (3.9), the quantum Potts 
Hamiltonian takes the following form: 



L-l q-1 



W = -Ef E^U 



1=1 



a 

1=1 H p 



Dual operators which satisfy the Potts algebra in 
Eq. (3.4) are defined as follows: 

Mi = fii_iftJ (l = 2,L), 

Mi = nj 

L-2-1 

ni= n M w M i+i+x {1 = 1, L-i), 

i=0 

n L = M L . (3.11) 

Under the duality transformation the Hamiltionain in 
Eq. (3.10) is changed into 



9-1 i-i r 3-1 



9-1 



"—E^E^-E^B^r-vE^- 

i=2 ^ p=l i=l 1 p=l q p=l 

(3.12) 



Having fixed boundary conditions with Sl = 1 amounts 
to have a vanishing transverse field T^, which eliminates 
the last anomalous term in the dual Hamiltonian. The 
coupling Ji-i is transformed into a transverse field act- 
ing on the Potts state at I while the transverse field T[ 
becomes a coupling between Potts states at / and 1 + 1. It 
follows that the original system with free-boundary con- 
ditions at I = 1 and fixed-boundary conditions at I = L is 
transformed into a dual system with fixed-boundary con- 
ditions at I = 1 and free-boundary conditions at I = L. 
The Hamiltonian is strictly self-dual when T; = J; either 
for homogeneous couplings T = J or, with inhomoge- 
neous couplings, when the sequence is symmetric. Then 
the self-duality point Ji = Ti corresponds to the critical 
point. 

In order to study the effect of a homoge- 
neous longitudinal magnetic field, the Hamiltonian in 
Eq. (3.10) can be completed by the diagonal field term 

-(fcA/9)£f=~iV(«J-l) -!]/(?- !)■ 

In the calculations we consider the behavior of the sys- 
tem in the vicinity of the self-duality point, taking the 
parametrization 



Ji = Ti 1 



t 

71 



(3.13) 



where t plays the same role as the temperature in the 
classical model: for t > (t < 0) we are in the bulk para- 
magnetic (ferromagnetic) phase with a vanishing (non- 
vanishing) order parameter. The local order parameter, 
i.e., the magnetization profile at site I, is defined as 
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mi 



q(0\5( Sl - 1)\0) - 1 



(3.14) 



where |0) is the ground-state eigenvector of H. On the 
other hand, the free-energy density of the 2D classical 



model is proportional to the ground-state energy per site, 
eo = (0\H\0) 22 

In the large-g limit one may use a systematic expan- 
sion in powers of l/^/q- To leading order the following 
eigenfunctions are degenerate: 



IV>4> 



|11. ..11), 

-L(|21...11> + |31...11> + . 



1 



|321 ... 11) + |421 ... 11) 



. + | l...ll)), 

+ |o(o-l)l...ll)], 



[|4321 ... 11) + |5321 ... 11) + ... + \q(q - l)(o - 2)1 . . . 11)], 



IV>l) 



1 



7 (L-l)/2 



\L . . .4321) + . . . + \q(q - l)(q - L) . . . (q - L + 2)1)], 



(3.15) 



all having the same diagonal contribution (ipi\H\il>i 



J2l=i T;, where we made use of Eq. (3.13). 



This degeneracy is lifted to the next order in 1/ ^fq and the secular problem is governed by the operator 

\ 



/ o ri 

r 1 tr 1 - h 



n 



i 



r 2 EUm-h) r 3 

r 3 EUm-h) 



V o 



(3.16) 



up to a constant term — J2f=i (L— l)h/ ^/q. One may 
notice that for the homogeneous problem with T; = T, 
the temperature and longitudinal field play the same role 
since we have the symmetry tY <-> —ft. The ground-state 
eigenvector of H is given by the linear combination 



|0> = J>hfc>, 



(3.17) 



i=l 



where Vi are the components of the leading eigenvector of 
-H in Eq. (3.16). Making use of Eqs. (3.14) and (3.15) 
we obtain the magnetization profile as 



mi 



(3.18) 



On the other hand, the parallel correlation length £m 
of the 2D classical model is related to the imaginary 
time spin-spin correlation function G;(r) of the quan- 
tum model. The long-time behavior of G;(r) has the 
asymptotic form 

G,(t) - (0|5[*,(0) - si(t)]\0) ~ exp(-rA) , (3.20) 

where A is the gap at the top of the spectrum of 
—H. Since A <~ one obtains information about 

the behavior of the correlations in the 2D classical sys- 
tem from the scaling properties of the spectrum of the 
Hamiltonian. 22 



in the large-g limit. In particular the surface magnetiza- 
tion is given by mi = \vi\ 2 and the bulk magnetization 
is defined as mj = m L / 2 - 

The perpendicular correlation length £_l gives the 
width of the interface region and, in the ordered phase, 
it can be defined through the relation 



mi\ 



(3.19) 



IV. SURFACE-INDUCED DISORDER IN 
REGULAR AND QUASIPERIODIC LATTICES 

A. Solution on the regular lattice 



For homogeneous interactions, i.e., when Ti = T, the 
complete solution of the Hamiltonian in Eq. (3.16) has 
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been given in Ref. 11. The critical exponents associated 
with the correlations parallel and perpendicular to the 
surface are given by 



v±. = v±_ = 



1 

3 ' 



-sf 



-100 



v\\ = v« = 



2 . (4.1) 




-50 50 

T z+l , T z+l 
tL , hL 



100 



I 1 /2ttZ 



(4.5) 




FIG. 1. Scaling plot of the singular part of the bulk energy 
density of the quantum Potts chain in the large-g limit for the 
regular (r = 1) and the Fibonacci lattice (r = 2, 3). The lines 
for finite systems with sizes L — 89, 233, 610, and 1597 are 
practically indistinguishable. 



For the surface magnetization one obtains 
x\ = 3 , 0i = l . 



(4.2) 



First, one may notice that conventional scaling relations 
in Eq. (2.12) and below Eq. (2.3) are satisfied. More 
importantly, the relations in Eq. (2.7), which are a con- 
sequence of the discontinuous behavior in the bulk, are 
also satisfied. 

Let us now consider the scaling of the free-energy den- 
sity in Eq. (2.4). For the quantum chain it corresponds 
to the singular part of the ground-state energy density, 
ep lng (i, h, L) — eo(t, h, L) — e (0, 0, oo), with the following 



finite-size scaling behavior: 25 

; (t, h, L) = L- {z+ ^e {tL z+1 , hL z+1 ) 



(4.3) 



According to analytical results [see Eq. (26) in Ref. 11] 
this relation is indeed satisfied for the homogeneous sys- 
tem with z — 2. The numerical results for the scaling 
functions are shown in Fig. 1. They agree with the result 
of a perturbation expansion eo(y, 0) = n 2 — y/2 + 0(y 3 ) 
to linear order. The scaling form of the magnetization 
profile in Eq. (2.5) is given by 



mi(t,h,L) =m l/L (tL z+i , hL 



-z+l\ 



(4.4) 



It is presented in Fig. 2 for h = and t < 0, t > 0, and 
t — 0. In the latter case the asymptotic scaling form is 
given by 



FIG. 2. Magnetization profiles of the homogeneous system 
with size L — 144 at different temperatures for free- and 
fixed-boundary conditions. The critical profile (t = 0) is the 
analytical result of Eq. (4.5). The interface for t < (t > 0) 
is localized at the free (fixed) boundary. 

In the ordered regime one can identify the interface re- 
gion with width £j_ ~ t _1 /( 1 + z ) near the surface at I = 1. 
In the high-temperature phase the interface has a width 
It is localized near the ordered surface at I = L. The 
two widths are related through = t), t > 0. 



2.5 



2.0 



O 1.5 



l.O 



0.5 



0.0 









V /v J\ A V =3 




Y a V =2 i 




/ V =l \ \a 



-500 



-250 



tL 



P +1 



250 



500 



FIG. 3. Scaling plot of the specific heat of the quantum 
Potts chain in the large-g limit for the regular (r — 1) and the 
Fibonacci lattice (r = 2,3). The lines are for finite systems 
with sizes L = 89, 233, 610, and 1597. For the Fibonacci 
lattice a double-peak structure is obtained. 
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According to Eq. (4.3), the specific heat of the system, 
c v ~ d 2 eo/dt 2 , diverges with L as c v <~ L 3 at the transi- 
tion point, a behavior which is compatible with Eq. (2.8) 
with z = 2. Outside the transition point the scaling 
function has been calculated numerically and is shown in 
Fig. 3. 



150 




-100 -50 50 100 

tL z+ \ hL z+l 



FIG. 4. Scaling plot of the energy gap of the quantum Potts 
chain in the large-q limit for the regular (r = 1) and the Fi- 
bonacci lattice (r = 2,3). The lines for finite systems with 
sizes L — 89, 233, 610, and 1597 indicate a good data col- 
lapse. The inset shows log-periodic oscillations of the scaled 
gap at larger values of the scaling variable for r = 2 and 
L = 89. 

Finally, we checked the scaling relation for the parallel 
correlation length in Eq. (2.9). It is related to the scaling 
behavior of the gap of the quantum Potts chain which is 
given by 



A(t,h,L) = L- z A(tL z+ \hL 



(4.6) 



The results of numerical calculations for the scaling func- 
tion A(j/,0) = A(0, — y) are shown in Fig. 4. They agree 
with the analytical result A(0,0) = 37r 2 [sec Eq. (21) in 
Ref. 11]. 



B. Fibonacci lattice 

The ID Fibonacci lattice is composed of two units A 
and B, and starting from A the sequence is generated 
through substitutions A — » AB and B — > A. The cou- 
plings or transverse fields on the lattice are also two val- 
ued, and Tb, depending on the letter in the sequence 
at the given position. The strength of the quasiperiodic 
perturbation is measured by the ratio r = Tb/^a- It 
has been noticed in Refs. 26 and 27, that the 2D layered 
Fibonacci lattice and the 2D Penrose quasilattice have a 
similar structure. On the Penrose lattice, starting from a 
straightlike surface, two different types of parallel layers 
can be identified, which arc distributed according to the 



Fibonacci sequence. Consequently the SID phenomenon 
studied here on the layered Fibonacci lattice should be- 
have in the same way on the Penrose lattice. 

Considering the Potts model in the largc-g limit 28 
the effective Hamiltonian in Eq. (3.16) at the transi- 
tion point t = h = takes the form of a ID quasiperi- 
odic Schrodingcr Hamiltonian, the properties of which 
have been extensively studied. 29 In particular, the scal- 
ing properties of the spectrum of Ti at special points, such 
as at the edge of the spectrum, are exactly known. 30 




FIG. 5. Finite-size-scaling exponents of the specific heat 
x Cv and surface magnetization X i . clS Qu function of the ampli- 
tude ratio r. The solid line represents the exact scaling result 
x Cv = z + 1 with z given by Eqs. (4.8) and (4.9). The dotted 
line is a guide to the eye. 

They are relevant for the strongly anisotropic scaling 
in our problem. The integrated density of states D(AE), 
measured from the top of the spectrum of — H has the 
scaling property 30 



D(AE) ~ \AE\ a , AE^O 



Here 



1 

a 



In A 



In 



A2 ' 



(4.7) 



(4.8) 



where <fi — (1 + V5)/2 is the golden-mean ratio, A = 
{8J-l+[(8 J-l) 2 -4] 1 /2}/ 2 with J = [3+(25+16/) 1 / 2 ]/8 
and I = (r + r _1 ) 2 /4. 

A relation between the exponent a and the anisotropy 
exponent z can be obtained by noticing the scaling rela- 
tion for the integrated density of states D' = D(AE') = 
bD(AE) when lengths are rescaled as V = L/b, since 
D(AE) is the number of states per unit length. Using 
also the scaling relation for the gaps AE' — b z AE, as 
given in Eq. (4.6), we arrive at the result 



1 

z = — 
a 



(4.9) 
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Consequently, the anisotropy exponent of the Potts 
model on the Fibonacci lattice is exactly known. Thus 
we can check the bulk scaling relations derived for SID, 
which contains only the anisotropy exponent z as a pa- 
rameter. The scaling behavior of the energy density in 
Eq. (4.3) with varying temperature and field is illustrated 
in Fig. 1. An excellent data collapse is obtained for the 
lengths indicated, i.e., for Fibonacci numbers F n with 
an odd index n. A different scaling function is obtained 
for n even. The scaling functions corresponding to the 
temperature and field perturbations are no longer re- 
lated through inversion, as they were for the homoge- 
neous problem. 




-30 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

-25 -20 -15 -10 -5 

ln\t\ 

FIG. 6. Temperature dependence of the gap (A) and the 
surface magnetization (mi) of the Potts model on the Fi- 
bonacci lattice with L — 1597 for different values of r, on a 
double-logarithmic scale. The numerically calculated points, 
which show log-periodic oscillations, are compared to asymp- 
totic scaling results shown by straight lines. The slope of the 
lines are v« = z/(z + 1) and f3\ = xi/(z + 1), for A and mi, 
respectively, where xi is taken from Fig. 5. Deviations from 
the linear behavior at small \t\ values are due to Unit-size ef- 
fects. Some sets of data were shifted vertically in order to 
improve the legibility. 

We have also studied the finite-size scaling behavior of 
the specific heat at the transition point which is given by 
Cy(t = 0, h = 0, L) ~ L Xa v and, according to Eq. (2.8), 
x Cnj = z+l. The exponents x Cv (r) were obtained through 
the extrapolation of two-point approximants using the 
Bulirsch-Stoer (BS) algorithm. 31 ' 32 As shown in Fig. 5 
this scaling relation is indeed satisfied in the whole range 
of amplitude ratio r investigated. 

The scaling form of the specific heat in Eq. (2.8) is 
shown in Fig. 3. For the Fibonacci lattice the scaling 
function around t = has a double peak structure with 
a distance between the peaks increasing with r. 

The finite-size scaling behavior of the energy gap given 
in Eq. (4.6) is tested numerically in Fig. 4 as a function 
of temperature and field. The unexpected maximum of 
the scaled gap at the transition point can be attributed 



to the presence of log-periodic oscillations, as shown in 
the inset for larger values of the scaling variable tL z+1 . 
Outside the transition point, the gap remains finite with 
an asymptotic temperature dependence A(t) ~ , with 
u\\ = z/(z + 1). The numerical results for the tempera- 
ture dependence, which is perturbed by log-periodic os- 
cillations, is shown in Fig. 6 for different values of r. 

Next we present the results about the surface scaling 
properties of the Potts model on the Fibonacci lattice. 
The determination of the order parameter, even at the 
surface, requires knowledge of the ground-state eigenvec- 
tor of H as shown in Eq. (3.18). With lack of an exact 
knowledge of this eigenvector we calculated numerically 
the surface magnetization at the transition point as a 
function of the size of the system, and the scaling expo- 
nent x\, defined in Eq. (2.11), was obtained through the 
extrapolation of two-point approximants using the BS al- 
gorithm. It is shown in Fig. 5 as a function of r. The 
surface magnetization was also calculated as a function 
of t for L = 1597 in order to check the scaling relation in 
Eq. (2.12) for the exponent [3\. The results are shown in 
Fig. 6. 




l/L 

FIG. 7. Magnetization profile of the quantum Potts chain 
at the transition point calculated on finite Fibonacci lattices 
with fixed- and free-boundary conditions. In the insets the 
first and second finite differences are taken for L — 2584. 

We close this section by presenting in Fig. 7 the mag- 
netization profiles calculated at the transition point on 
finite Fibonacci lattices. At first sight the shape of 
the curves for different lengths looks similar; however, 
a closer inspection shows that with increasing size there 
are more and more points where a stcplikc behavior de- 
velops. To look for a possible nonanalytical behavior of 
the curve in the thermodynamic limit we have calculated 
the finite size differences Aim; = L(m; + i — m;_i)/2 and 
A.2TTI1 = L 2 (mi + i + mi-i — 2mi). As seen in the insets 
of Fig. 7, the derivatives show a self-similar, multifrac- 
tal structure. With increasing size L the heights of the 
peaks are increasing, indicating divergent first and sec- 
ond derivatives. Consequently mi/L is an example of 
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a physical observable with nonanalytical mathematical 
properties. 



V. DISCUSSION 

In this paper we have presented a scaling picture for 
bulk first-order transitions associated with SID, i.e., with 
a continuous surface transition. We have shown that fluc- 
tuations of the interface separating the ordered (bulk) 
and disordered (surface) regions govern the scaling be- 
havior of bulk quantities. Power-law singularities are 
present and the bulk free-energy density shows a strongly 
anisotropic scaling behavior at the first-order transition 
point, which can be associated with an anisotropic DFP. 

The discontinuous nature of the bulk transition fixes 
the form of the scaling exponents, which can all be ex- 
pressed in terms of the anisotropy exponent z and the 
dimension of the system. Also the functional form of the 
scaling functions is particular. For example, response 
functions at the transition point have a power-law-type 
size-dependent peak, from which a (5-function singularity 
develops in the thermodynamic limit. 

The scaling picture for SID has been tested on the 
2D Potts model in the large-g limit, both on the regu- 
lar lattice and on the layered Fibonacci lattice (the lat- 
ter is related to the 2D Penrose quasilattice). For the 
Fibonacci lattice the anisotropy exponent z varies con- 
tinuously with the amplitude ratio of the apcriodicity 
in an exactly known way. The study of the associated 
variation of other scaling exponents provides us with a 
stronger test of the SID scaling relations. 

The marginal nature of the quasiperiodic perturbation 
of the Fibonacci type should be explainable through a 
relevance-irrelevance criterion, like the Harris criterion 13 
for critical systems in the presence of quenched disorder. 

For second-order transitions the relevance or irrele- 
vance of aperiodic perturbations is generally related to 
the size dependence of the fluctuations of the sequence, 
which is measured by the wandering exponent lo. For a 
ID variation, as considered in this paper, lo is defined 
through the cumulated deviation from the average 33 



A r (x n ) = £(r* - r) 
1=1 



(5.1) 



where T = limi, n _ >(X) (l/L„) J2i=i is the average field 
and L n is the length of the sequence after n substi- 
tutions (i.e., the Fibonacci number F n ). As shown in 
Refs. 34 and 35, for a thermal aperiodic perturbation, 
the crossover exponent at the fixed point of the pure sys- 
tem is ip = 1 + v±(io — 1), where v± is the perpendicular 
correlation length exponent of the pure system. Such a 
perturbation is relevant (irrelevant) when ip > (ip < 0). 

A naive application of this criterion to our problem, 
with v±_ = 1/3 according to Eq. (4.1), leads to a marginal 



value of the wandering exponent lo = — 2. Thus the per- 
turbation should be relevant for the Fibonacci sequence 
with lo = — 1. 

Actually, the self-dual form of the aperiodic perturba- 
tion considered in this work, with J/ = T; in Eq. (3.10), is 
weaker than a standard thermal perturbation for which, 
for example, Ji fluctuates aperiodically while T; is kept 
constant. 

In order to obtain a relevance-irrelevance criterion 
adapted to our problem, let us consider the form of the 
Hamiltonian matrix in Eq. (3.16) and compare the per- 
turbation associated with a deviation t from the transi- 
tion point, on the one hand, to the perturbation intro- 
duced by the aperiodic fluctuations on the other hand. 

The temperature perturbation contributes to the diag- 
onal elements by terms of the form 



(5.2) 



According to Eq. (5.1) when k is large the last term be- 
haves as tk u and may be neglected since the wandering 
exponent is always smaller than 1. Taking the average 
at a length scale L > 1 of the leading contribution to 
Eq. (5.2), one obtains the following estimate for the tem- 
perature perturbation: 



Tt 



L-l 



5 t (L) = -J2k~tL. 



(5.3) 



fe=i 



The aperiodic fluctuations contribute mainly to the off- 
diagonal terms and using Eq. (5.1), their mean value at 
the length scale L is given by 



S r (L) = yA r (L) ~ L^ 1 



(5.4) 



The appropriate length scale is the perpendicular corre- 
lation length £j_ and the relevance of the aperiodic per- 
turbation depends on the value of the ratio 



(5.5) 



where we used Eq. (2.2). Thus for self-dual aperiodic 
perturbations the crossover exponent is 



<Psd 



l + v±(w-2) . 



(5.6) 



It follows that with v± = 1/3 such perturbations are 
marginal when lo = — 1, in agreement with our results 
for the Fibonacci lattice. 

Finally, let us briefly mention that we have also stud- 
ied the effect of an aperiodic perturbation following the 
period-doubling sequence. 36 This sequence is generated 
through the substitutions A — ► AB,B — > A A, and its 
wandering exponent is lo — 0. According to Eq. (5.6) 
this perturbation is relevant at the first-order transition 
point of the g-state Potts model. Indeed, using the results 
of an exact renormalization group method in Ref. 16, 
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06 



the anisotropy exponent is found to be formally infinite. 
More precisely the correlation lengths are related through 

ln£||~#, (5.7) 
where fi = 1/2 for the period-doubling sequence. 37 
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